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A única fonte luminosa artificial que reproduz as condições da natu- 
reza proporcionando uma luz difusa com 5.500º Kelvin de temperatura 
de côr, ideal portanto para todos os trabalhos de fotografia a cores, 


Foram realizados ensaios utilizando cartas padrão Kodak, de côr e 
de cinzentos, e filme Ektachrome Daylight, os quais revelaram não 
existirem diferenças sensíveis entre diapositivos obtidos no exte- 
rior às 12 horas em dia descoberto, e diapositivos obtidos em interior, 
sem mistura de outras fontes luminosas e nas mesmas condições de 
velocidade de obturador e abertura de diafragma, mas sob a acção de 
luz directa da iluminação fluorescente OPTIMA. 


Peça hoje mesmo uma demonstração sem compromisso ao represen- 
tante exclusivo para todo o território Português. 
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administração actualizada, a capa- 
cidade técnica dos seus engenhe:- 
ros e o profissionalismo com que 
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suas responsabilidades, são valo- 
res reais e positivos que, de ma- 
neira decisiva, contribuem para 
o engrandecimento da Standard 
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Associada da International Tele- 
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Eléctrica fabrica numerosos pro- 
dutos, desde os minúsculos tran- 
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serviço publico, que são exporta- 
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O PAPEL DOS LABORATÓRIOS NA ENGENHARIA CIVIL 


RESUMO 


Procura dar-se uma visão de conjunto sobre a contri- 
buição que os laboratórios hoje dão para a resolução dos 
problemas de engenharia civil, É referido o estudo 
dos materiais, a averiguação das acções a considerar 
no projecto, a definição das formas a dar às obras, 
o estudo dos processos de construção a seguir, e a veri- 
ficação do comportamento das construções. Aborda-se 
ainda o problema dos tipos de laboratórios em face 
da satisfação das necessidades da indústria da cons- 
trução. 


por MANUEL ROCHA** 


SYNOPSIS 


The paper gives an overall view of today's contribus 
tion of laboratories to the solving of civil engineering 
problems. It concerns the study of materials, the inves- 
tigation of the actions to be considered in the design, 
the definition of shapes for the constructions, the study 
of the construction techniques to be followed, and the 
observation of the behaviour of the constructions. The 
problem of the types of laboratories is also tackled with 
a view to meeting the needs of the industry of con- 
struction. 


Sinto-me muito honrado e sensibilizado por ter sido distinguido pelo convite para falar nesta 
cerimónia dedicada a marcar a inauguração do Raymond E. Davies Hall. Compreenderão a emoção 
de quem vem do distante Portugal tomar a palavra, em tal momento, no país que ainda há pouco 
levou a cabo o maior feito do Homem e numa Universidade que é o maior conglomerado de saber 
que jamais existiu e que tanto está contribuindo para modelar o futuro deste país. Acresce que desde 
que iniciámos o nosso Laboratório em Lisboa — já lá vão mais de 20 anos — temos procurado desenvolver 
estreita cooperação com o Departamento de Engenharia Civil, o qual ao longo deste período tem 
mantido posição do maior destaque, quer na investigação, quer no ensino da engenharia civil. E é 
oportuno referir que o primeiro professor do Departamento que nos visitou foi o Professor Raymondi 
Davies a quem todos os de nós que têm tido o privilégio de o conhecer dedicam a maior admiração 
e estima pela sua tão rica e inconfundível personalidade. 


1. Quanto ao tema da minha oração — O papel dos laboratórios na engenharia civil — começarei 
por referir que tenho tido a oportunidade de observar esse papel em países com os mais diversos 
graus de desenvolvimento. Tenho verificado que, tal como aconteceu na evolução histórica da 
engenharia, em que o primeiro laboratório foi criado para servir a engenharia civil — trata-se do 
«Laboratoire des Ponts et Chaussées», de Paris, instituído em 1851 e destinado ao estudo de materiais 
de construção — , ainda hoje nos territórios em vias de desenvolvimento são os laboratórios destinados 
a apoiar as realizações do engenheiro civil aqueles cuja necessidade se faz sentir em primeiro lugar. 


* Discurso de inauguração dos novos Laboratórios do Departamento de Engenharia Civil da Universidade 
da Califórnia, Berkeley. 
** Director, Laboratório Nacional de Engenharia Civil, Lisboa. 
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E tal é consequência do desenvolvimento de um território ter de começar pela criação de infra- 
-estruturas, isto é, por obras de engenharia civil, as quais, pela sua natureza, não podem ser importadas. 
Têm de ser executadas localmente, assim se criando condições para o desenvolvimento de laboratórios 
de engenharia civil. À medida que um território se desenvolve cada vez o projecto e a execução de 
obras tira maior proveito da acção dos laboratórios. Considero, no entanto, que mesmo nos países 
mais evoluídos não são em regra plenamente aproveitadas as tão variadas possibilidades, que adiante 
referirei, oferecidas pelos métodos experimentais, mesmo no respeitante ao aproveitamento das 
capacidades dos laboratórios já existentes. 

Vejamos então em que é que os laboratórios podem ser úteis à engenharia civil, à qual incumbe 
a produção e conservação de edifícios e obras públicas e, por vezes, mesmo a exploração das obras 
construidas. 

Os laboratórios dão a sua contribuição em relação aos materiais a usar, às acções exteriores a 
considerar, às formas a adoptar, aos próprios processos de construção a seguir, e compete-lhes ainda 
a tarefa de verificação do comportamento da obra, isto é, do produto final da indústria da 
construção. 

Ao falar da contribuição dos laboratórios desejo esclarecer que nela incluo a variada gama de 
actividades que vai desde o simples ensaio de rotina até à investigação destinada a servir a engenharia 
civil. E é oportuno frisar que cada vez melhor se consegue planear e conduzir a acção de investigação 
de modo a desta derivarem, mesmo a curto prazo, reais benefícios. Acresce que o trabalho de 
laboratório, isto é, a experimentação, conduz naturalmente à investigação pois que, como afirmou de 
maneira tão incisiva o Professor italiano Arturo Danusso, «... cabe ao experimentador a honra 
de receber directamente a voz da Natureza para depois a interpretar. Esta voz, quando escutada com 
humildade, ensina muitas coisas, rectifica previsões e destrói esquemas prévios, ao mesmo tempo 
que sugere outros». 


2. O estudo dos materiais, que, até fins do século passado, constituía a quase exclusiva 
atribuição dos laboratórios, continua sendo hoje uma das mais volumosas e frutuosas actividades dos 
laboratórios destinados a apoiar a engenharia civil. Na verdade, são os laboratórios que permitem 
caracterizar as variadíssimas propriedades dos materiais de que depende o sucesso das obras, tais 
como a resistência sob acções estáticas e dinâmicas, a permeabilidade a líquidos e gases, as variações 
de dimensões resultantes de variações de temperatura e humidade e de acções químicas, as propriedades 
térmicas, as propriedades acústicas, a resistência ao fogo, a durabilidade sob a acção de agentes 
químicos e biológicos, e mesmo as características estéticas, como a cor. 

Muitos ensaios, em especial os ensaios mecânicos, visam a determinação de propriedades a 
partir das quais, por meio de teorias, é possível prever o comportamento dos materiais em obra. 
Quando não se dispõe de teoria de previsão adequada, tem de se lançar mão a ensaios que simulem 
as acções a que os materiais se encontrarão submetidos em obra, quer as provenientes do meio físico 
ambiente, quer do utilizador. Por vezes tem mesmo de se recorrer a ensaios «in situ», em especial 
para ser possível ensaiar volumes suficientemente grandes para serem representativos, como acontece 
com frequência quando há a tratar problemas respeitantes a solos e rochas. 

No que respeita aos materiais, podem distinguir-se os três seguintes tipos de contribuições dos 
laboratórios: i) os ensaios destinados à prévia selecção e definição de características dos materiais 
a usar em cada obra, em particular dos materiais como o betão e materiais para aterros, que serão 
fabricados no estaleiro, isto é, sob a responsabilidade do engenheiro civil; ii) os ensaios de recepção 
de materiais durante a construção, quer dos produzidos em fábrica, quer em estaleiro ; iii) a investigação 
destinada a conhecer e a melhorar as propriedades dos materiais disponíveis e a descobrir novos 
materiais. 

Apesar de nas obras do engenheiro civil ainda persistir a utilização de numerosos materiais 
tradicionais, verifica-se a introdução a ritmo acelerado, e que certamente se acentuará no futuro, de 
novos materiais, sobretudo na construção de edifícios, o que confere cada vez mais relevância ao 
papel desempenhado pelos laboratórios. Na verdade, esta explosão de novos materiais não permite 
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solicitações de carácter aleatório sobre as estruturas, e a busca de um melhor conhecimento dos 
fenómenos de erosão. 

No que respeita à hidráulica marítima, os laboratórios prestam larga colaboração, em especial 
no estudo da agitação nas zonas a proteger com obras portuárias, da acção das vagas sobre as 
estruturas, em especial molhes e esporões, da dinâmica e da evolução dos fundos de estuários, 
incluído o complexo problema da sedimentação das vasas, e, finalmente, referirei o problema da 
conservação das praias que está assumindo cada vez mais importância. 

No estudo dos problemas marítimos a ondulação era, até há poucos anos, simulada nos modelos 
por meio de ondas monocromáticas. Nos últimos tempos tem-se procurado uma simulação melhor 
ajustada à realidade, introduzindo nos modelos uma ondulação constituída pela sobreposição de ondas 
com diversas amplitudes e frequências, tal como na Natureza. Os modelos também estão permitindo 
a determinação das acções mecânicas, de carácter aleatório, exercidas sobre as estruturas marítimas. 

No respeitante aos problemas de agitação, estão sendo feitos progressos no tratamento analítico 
dos fenómenos de difracção, refracção e reflexão que estão permitindo simplificar e poderão vir a 
dispensar a realização de ensaios. 

Na hidráulica fluvial os modelos prestam também relevantes serviços no estudo dos problemas 
de controle de cheias, de navegação, de estabilidade de leitos, e de transporte sólido, por arrastamento 
e em suspensão. Alguns dos problemas de escoamento com fundo fixo podem já ser abordados com 
vantagem por via analítica, quando as condições não são complexas. No domínio dos escoamentos 
com transporte sólido, os modelos, devidamente calibrados por comparação com os protótipos, 
constituem uma ferramenta indispensável. No entanto, nos últimos anos, graças ao progresso no 
tratamento dos processos estocásticos, possibilitado pelos computadores, foi retomado, com perspectivas 
de êxito, o estudo dos fenómenos do transporte sólido em termos probabilísticos. 

Não quero deixar de mencionar a contribuição dos laboratórios para a resolução de problemas 
de engenharia sanitária, como os postos pelas redes de abastecimento de águas e esgotos, pelos 
dispositivos sanitários em edifícios, e pela poluição de águas por afluentes domésticos e industriais. 

No que respeita à contribuição dos modelos para a definição das formas e dimensões a dar às 
construções, referimos em primeiro lugar os modelos estruturais e depois os modelos hidráulicos, 
O projectista pode ainda tirar proveito de modelos puramente funcionais, isto é, destinados a averiguar 
se dada função de uma obra é devidamente desempenhada. Tais modelos têm especial interesse nos 
estudos de edifícios, domínio no qual há a satisfazer numerosas condições além da resistência 
estrutural. Assim, são usados modelos no estudo da acústica de salas, da insolação e iluminação 
natural, da ventilação, etc. 

Os modelos de quaisquer dos variados tipos que temos vindo a referir podem ser utilizados, 
quer para a verificação de dada solução prevista pelo projectista e colheita de sujestões para o seu 
aperfeiçoamento, quer para a tarefa de investigação de novas formas a dar às construções, sob os 
múltiplos aspectos que interessam ao engenheiro, quer ainda para a verificação de teorias de previsão 
do comportamento das construções. 


5. Dentro da concepção de investigação integrada que adiante referiremos, verifica-se que os 
organismos de investigação começam, nos últimos anos, a ocupar-se também dos problemas levantados 
pela própria construção das obras. A investigação está cobrindo não só a fase de execução mas 
também os problemas de gestão levantados pela indústria da construção. Mais do que isso, em certos 
países a actividade de investigação está-se mesmo estendendo à economia aplicada à engenharia, 
tratando-se problemas como o planeamento, num dado território, das realizações a empreender de 
modo a optimizar a utilização dos meios disponíveis, a satisfação de necessidades por meio de 
equipamentos evolutivos, como por exemplo no caso de edifícios para habitação e de estradas, a 
obsolescência dos equipamentos, etc. 

No respeitante à gestão estão sendo tratados todos os problemas relativos à execução da 
construção, como o planeamento da obra, a sua programação, a elaboração de propostas, a contabilidade 
de custos para previsão e verificação do custo da construção, e a produtividade global e de cada um 
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dos trabalhos da obra. Na investigação dos problemas da gestão estão sendo dados os primeiros 
passos para o desenvolvimento de um sistema de gestão integrada. 

Quanto à fase de execução, a investigação está cobrindo o desenvolvimento de equipamento e 
os métodos e técnicas de execução, áreas estas nas quais se oferecem amplas perspectivas de 
sucesso. 


6. Finalmente, após ter sido referida a participação dos laboratórios no respeitante aos materiais, 
as acções sobre as obras, às formas destas e aos processos de construção, vejamos qual a sua 
contribuição para a observação do comportamento das construções. 

A observação consiste na colheita sistemática de informações sobre o comportamento estrutural 
e funcional, com vista: i) a verificar a justeza das previsões contidas no projecto da obra em 
observação, e consequentemente, julgar o mérito do património de conhecimentos em que o projectista 
se apoiou; e ii) a obter novos conhecimentos, a aproveitar em futuras realizações. 

Os projectistas, construtores, utilizadores e entidades responsáveis pela exploração das obras 
participam, em graus vários conforme os casos, na tarefa de observação mas esta exige muitas vezes 
pessoal e meios experimentais especializados que obrigam à intervenção dos laboratórios, em particular 
quando a observação visa o alargamento de conhecimentos. 

Quanto à tarefa de observação, referiremos em primeiro lugar a observação do comportamento 
dos materiais, em especial a averiguação da evolução de propriedades no tempo que possa ter 
repercussões desfavoráveis no comportamento estrutural ou funcional. 

Uma das mais relevantes contribuições dos laboratórios para a observação das obras tem sido 
a observação do comportamento estrutural, a qual interessa por vezes empreender mesmo no decorrer 
da construção. De um modo geral, a apreciação do comportamento estrutural implica a observação 
dos efeitos das solicitações, isto é, deslocamentos, deformações e tensões, assim como a observação 
das próprias solicitações actuantes na estrutura tais como o peso próprio, pressões exercidas pelo 
vento, aceleração de sismos, assentamento de fundações e variações da temperatura. Além disso, 
a interpretação dos resultados da observação exige, em regra, a determinação, em laboratório, de 
diversas propriedades dos materiais da estrutura, em especial das propriedades mecânicas e termo- 
-higrométricas. 

Como apreciação global da actividade de observação do comportamento estrutural conduzida 
pelos laboratórios, considero que na maior parte dos domínios tem sido insuficiente. É particularmente 
magra a informação de que se dispõe sobre o efeito de solicitações cuja probabilidade de ocorrência 
é baixa, tais como sismos e ventos de elevada velocidade. Penso que a observação do comportamento 
estrutural devia ser empreendida mais largamente, tanto mais que a tarefa de observação está 
presentemente muito facilitada pelos desenvolvimentos realizados nos últimos anos no domínio da 
aparelhagem de medida. A intensificação da observação concorrerá não só para o progresso dos 
conhecimentos mas também para o amadurecimento da formação dos intervenientes no projecto e 
execução das estruturas, dada, diremos, a visualização que proporciona dos fenómenos do funciona- 
mento estrutural. 

Terminarei a referência que tenho estado a fazer à contribuição dos laboratórios para a 
observação das construções mencionando ainda a observação do comportamento funcional, isto é, 
respeitante à utilização da obra. Neste domínio há numerosas observações que podem ser feitas 
pelos utilizadores e pelas entidades que exploram as obras, tais como a observação da adequação 
dos espaços de um edifício às funções que este deve desempenhar, das infiltrações numa obra 
hidráulica ou do tráfego numa estrada, mas há muitas outras observações para cuja realização os 
laboratórios estão especialmente preparados. Será, por exemplo, o caso da observação quantitativa 
das condições de conforto oferecidas por um edifício, da derrapância e outras características do 
pavimento de estradas, etc. 

Acabámos de passar em revista as contribuições dos laboratórios para a solução dos problemas 
da engenharia civil. Apesar de não ter podido senão abordar superficialmente o tema, acho que ficou 
bem evidente a extrema diversidade dessas contribuições e a sua relevância. 
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7. Será oportuno referir brevemente os tipos de laboratórios que prestam apoio à engenharia 
civil e o papel que desempenham. 

Menciono em primeiro lugar os laboratórios que se ocupam sômente de ensaios de rotina. 
Entre estes há a considerar os laboratórios comerciais de ensaios de materiais e de modelos e os 
laboratórios de estaleiro destinados a apoiar a fiscalização da execução das obras. 

O outro tipo de laboratório, que desejaria analisar com certo pormenor, é aquele em que é 
empreendida actividade de investigação, associada por vezes a actividade de desenvolvimento. São 
muito diversas as índoles dos laboratórios de investigação no respeitante a política de investigação, 
extensão do domínio coberto, dimensão e vinculação. 

No tocante à política de investigação seguida, classificarei os laboratórios em dois grupos: os 
que se ocupam de dadas disciplinas e os que empreendem investigação que designo por integrada. 

A investigação conduzida nas Universidades é em regra de natureza disciplinar, o que resulta 
de facto de as Universidades terem, naturalmente, como células da sua estrutura as disciplinas do 
conhecimento. 

A investigação integrada, isto é, aquela que visa a solução de problemas com apoio em todos 
os conhecimentos disponíveis, independentemente das disciplinas envolvidas, é em regra conduzida 
por organismos do sector público, pois na indústria da construção não se verificam condições 
favoráveis à existência de órgãos de investigação integrada no sector privado, contrariamente ao que 
ocorre na maior parte das restantes indústrias. Dada a grande extensão e diversidade dos campos 
cobertos pela engenharia civil, a investigação integrada tem de ser conduzida nos países com maior 
dimensão económica em laboratórios que cobrem sômente um ramo da indústria da construção, tal 
como edifícios ou estradas ou aproveitamentos hidráulicos. No entanto, em países com menor 
dimensão económica, como Portugal, considero possível concentrar a investigação na engenharia 
civil, conduzida de maneira integrada, num só laboratório de grande dimensão, diremos como um 
quadro da ordem de um milhar de pessoas. 

Dentro do conceito de investigação integrada, os organismos que dela se ocupam tanto 
empreendem trabalho de natureza experimental como teórica e por isso será preferível designá-los 
por institutos em vez de laboratórios, apesar de predominar, no caso da engenharia civil, a actividade 
de natureza experimental. 

Estou considerando esta matéria dos tipos de organismos de investigação porque ela é de 
grande importância do ponto de vista do tema que estou tratando, isto é, do papel dos laboratórios 
na engenharia civil. Na verdade, um instituto de investigação integrada, pela sua própria natureza, 
tem de manter actualizada a apreciação global das necessidades presentes e futuras do ramo da 
indústria da construção que é seu mister apoiar, tem de procurar dispor dos meios humanos e 
materiais exigidos pela sua satisfação e não pode deixar de se preocupar com a efectiva utilização 
dos resultados da sua actividade, quer através de acções de difusão, quer de desenvolvimento. 
A referida apreciação global de necessidades tem o maior interesse pois permite uma repartição judiciosa 
da capacidade de investigação disponível, evitando que seja dada atenção excessiva a certos problemas 
e, sobretudo, a omissão de outros que podem não ser atraentes mas que são importantes. 

A investigação integrada facilita o estabelecimento de uma associação íntima das actividades 
de investigação com a prática, sobretudo se os institutos em que é empreendida, além da actividade 
de investigação planeada, também se prestem a dar a sua colaboração na solução de problemas postos 
individualmente por projectistas e construtores, e mesmo na realização de ensaios de rotina. Penso 
ser esta a boa solução, desde que sejam respeitadas certas proporções, que variam de caso para caso, 
entre os volumes dos vários tipos de actividades. 

O diálogo permanente entre a investigação e a prática tem duas consequências com a maior 
projecção. Em primeiro lugar, ele é indispensável para os organismos de investigação poderem 
adquirir uma compreensão clara do contexto técnico, económico e social onde se deve inserir a sua 
acção, condição básica para a eficiência da investigação, desde a formulação dos problemas a resolver, 
à condução dos trabalhos, ao estabelecimento de conclusões e às medidas visando assegurar a 
utilização dos resultados. Em segundo lugar, o referido diálogo contribui para uma mentalização dos 
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projectistas e construtores que os leva a poderem dar uma melhor contribuição para o estabelecimento 
de planos gerais de investigação e a tomarem uma posição mais receptiva perante os resultados da 
investigação. Além disso, essa mentalização contribuirá para evitar que sejam desperdiçadas as 
oportunidades que surgem no projecto e execução de cada obra para os institutos de investigação 
darem a sua colaboração na solução de problemas. 

Dados os méritos da investigação integrada, considero que devem existir, em cada país, 
institutos de investigação integrada destinados a servir os vários ramos da indústria da construção 
e que, além disso, se deve procurar intensificar o grau de integração da investigação conduzida nas 
universidades. É no domínio dos edifícios, em especial dos destinados a habitação, que no momento 
presente mais se faz sentir em muitos países a necessidade de ser dado mais desenvolvimento à 
investigação integrada. 


8. Ao falar-se em investigação põe-se o problema do grau de desenvolvimento que é pertinente 
dar-lhe. Não há muitos decénios tal questão não era formulada em termos globais, mesmo nos países 
mais avançados, deixando-se a sua evolução às iniciativas individuais. Mas, presentemente, dentro 
do crescente esforço da Humanidade para se determinar a si mesma, verifica-se um grande esforço 
para definir critérios que permitam estabelecer tão objectivamente quanto possível, dentro da óptica 
custo-benefício, o investimento que deva ser feito na investigação. Trata-se, como é óbvio, de 
problema muito complexo que implica a consideração, em cada país, de directrizes políticas, de planos 
de acção económica, da disponibilidade de quadros e ainda de outros parâmetros. 

No caso da engenharia civil, verifica-se em países da Europa Ocidental, como a Grã-Bretanha 
e Portugal, que a despesa global com a investigação, do sector público e privado, atinge cerca de 0,5" 
do valor total das obras executadas, isto é, edifícios, incluídos os de habitação, e obras públicas. 
Nos Estados Unidos o valor total das obras executadas atingiu, em 1968, cerca de 85 000 milhões de 
dólares mas não consegui obter o valor da despesa global com a investigação. Se se admitir que esta 
atinge também 0,5º/0 do valor das obras executadas, tal corresponderia a uma despesa anual com a 
investigação de 425 milhões de dólares. Este valor corresponderia a cerca de 40 000 pessoas envolvidas 
na investigação, supostas em «full-time». Suponho que este número não é atingido. 

Os números apresentados mostram a importância que atinge hoje no Mundo o valor dos 
investimentos na investigação que serve a construção e portanto a necessidade de cuidar da eficiência 
da acção de investigação, não só pelo montante da despesa mas sobretudo pela grande repercussão 
que essa acção pode ter na indústria da construção. A investigação deve ser encarada tal como uma 
actividade industrial que tem de escolher cuidadosamente os produtos a fabricar, os métodos de 
produção a seguir e que não pode deixar de cuidar da venda dos seus produtos. 

Quanto à fixação do grau de desenvolvimento a dar às actividades de investigação, creio que 
será hoje possível, graças à volumosa informação disponível, estabelecer para cada um dos grandes 
ramos da indústria da construção o investimento que convirá fazer, o qual variará de país para país. 

Penso que na maior parte dos países se justifica, dentro de uma concepção custo-benefício, 
incrementar a actividade dos laboratórios que apoiam a indústria da construção, uma vez que a 
investigação seja conduzida aproveitando todo o vasto património de conhecimentos e experiência 
de que hoje se dispõe para assegurar a sua eficiência, em especial dentro da concepção de investigação 
integrada. A necessidade de expansão dos laboratórios é consequência da profunda modificação que 
se está operando na engenharia civil, a qual está passando rápidamente de uma actividade essencialmente 
empírica, como se verificava ainda há alguns decénios, a uma actividade comandada pelo conhecimento 
científico. 
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NOTE PRELIMINAIRE SUR 


RESUMO 


Apresenta-se o processo geral que condusirá à defi- 
nição de um indice de alterabilidade a partir de ensaios 
de alteração acelerada realizados em laboratório sobre 
os minerais essenciais das rochas igneas e também sobre 
os principais tipos destas rochas com granularidades 
e texturas diferentes, 


On a défini un indice de qualité micropétro- 
graphique (MENDES, AIRES-BARROS & RO- 
DRIGUES, 1966), qui permet de: 


a) caractériser I'état d'altération des massifs 
rocheux ; 

b) obtenir la corrélation entre cet état géo- 
technique et quelques grandeurs mécaniques 
(E et 7) des roches de ces massifs rocheux dont 
on essaie des échantillons au laboratoire ; 

c) expliquer les valeurs de E et 7 d'accord avec 
la composition minéralogique et la texture des 
roches que Ion étudie au microscope polarisant. 


Tout en se basant sur I[analyse modale 
on a défini k, lindice de qualité micropétro- 
graphique par : 


n 
>, P; Xi; 


PARES de À ER 


E ã 
o P; Y; 


j= 


ER 


dont les n valeurs de x, sont les pourcentages 
modales des minéraux «sains», avec incidence 
favorable dans le comportement mécanique des 


UN 
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* 
INDICE D'ALTERABILITE 
par L. AIRES BARROS ** 
Laboratoire de Minéralogie et Pétro- 
logie, Institut Supérieur Techniques 
Lisbonne-Portugal 
SYNOPSIS 


The general process for definition of an alterability 
index is presented, 

This index will be obtained after several laboratory 
aging tests on essential minerais of igneous rocks and 
also on the principal types of these rocks with different 
textures and granularities. 


roches, et dont les m valeurs de Y, sont les 
pourcentages modales des minéraux altérés ou 
des minéraux phyliteux qui montrent disposition 
orientée dans les roches, et aussi des microfissures 
et des vides, tous avec incidence défavorable 
dans le comportement mécanique des roches. 

Tout récemment indice k a été appliqué à 
étude de l'anisotropie des roches endogênes. On 
a obtenu de bonnes corrélations entre les élli- 
psoides des indices micropétrographiques et les 
éllipscides des tensions de rupture et les élli- 
pscides des modules d'élasticité (RobRiIGUES & 
AlIRES-BARROS, 1970). 

Le degré d'altération d'une roche (m) est 
fonction de: 


me) .. cus (2) 


ou: 


i==f, (n,s) et e = f, (9, ph, Eh, 1,b) 

i, représente les facteurs intrinsegues à la roche, 
comme le type des minéraux (n) et la surface 
exposée à lataque (s) qui est fonction de la 
texture et de la microfissuration et microfractu- 
ration de la roche et aussi de la quantité de vides. 
e, représente les facteurs extrinseques, comme la 


* Comunicação apresentada na IV secção (Weathering phenomena of soils and rocks) do 1.º «International Con- 
gress of the International Association of Engineering Geology», realizado em Paris, sob os auspícios da UNESCO, 


de 8 a 11 de Setembro último. 


** Trabalho realizado dentro do âmbito do projecto de investigação «Geoquímica da alteração experimental 
de rochas», subsidiado pelo Instituto de Alta Cultura do Ministério da Educação Nacional. 
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température 1, le ph, le Eh, le débit | des fluides 
de percolation et encore les forces biotiques b. 

L'altérabilité (m,) sera fonction d'un «temps 
humain» (t) qui est fondamental dans l'étude 
que nous concerne, vu que le temps de I'altération 
est principallement un «temps géologique» (T). 


Ainsi: 
m =F(iet,T)..... (3) 
Si on considere que: 


a) k=F (n,s,t,T); 

b) les facteurs extrinsêgques e, peuvent être 
postulés au laboratoire par des essais accélerês 
d'altération des roches (essais de percolation 
libre intermitente, essais de vieillissement par 
agitation intermitente, essais de percolation 
forcée, essais de percolation dans des extracteurs 
de Soxhlet, essais de vieillissement par irradiation 
thermique, etc); on pourra chercher un indice 
d'altérabilité défine par: 


d k , 4 
k' = — significativement different de zéro 
dt 


(4) 


Les essais de laboratoire en cours, en ce mo- 
ment, fixent les facteurs extrinseéques e pour: 

a) les minéraux éssentiels des roches ignées, 
voire les components des séries de Goldich 
(1938): 

b) les principaux types de roches ignées (gra- 
nites alcalins, granites calcoalcalins, syénites, 
syénites néphéliniques, diorites, gabbros, basaltes 
à olivine et dolérites); 

c) les trois principaux types de granularité 
(grain grossier, moyen et fin) d'un même type de 
roche ignée (syénite néphelinique). 

A partir d'un grand nombre d'essais d'altéra- 
tion accélerée au laboratoire, dont on a fait 
brieve mention, on peut mesurer les vélocités 
rélatives d'altération (pour des conditions exo- 
génes définies à Iavance), des minéraux des 
séries de Bowen. 

De même on étudie les taxes d'altération arti- 
ficielle des principaux types des roches ignées 
(déjá mentionnées) et on determine l'influence de 
la grannularité, soit de la surface espécifique à 
ataque par les agents exogênes. 


10 


Les velocités rélatives d'altération seront deter- 
minées au moyen des essais physiques (varia- 
tions du pH et de la résistivité) et chimiques 
(enrichissement en E Nel. CT Re, Mg?” 
et en silice) des fluides de percolation et aussi des 
essais physiques (variations de porosité, de per- 
méabilité à lair et de susceptibilité magnétique, 
pertes de poids), chimiques (variations de com- 
position) et minéralogiques des minéraux et des 
roches en étude. 

Ces essais doivent permettre de determinar les 
incidences de la composition minéralogique (n) 
et de la granularité et microdiscontinuités (s) 
dans la velocité d'altération. Les expériences 
accomplies sont encore en nombre insuffisant 
pour présenter des conclusions déja définiti- 
ves. 

Mais supposons que, pour la série de miné- 
raux em étude, on obtiendra la séguence de 
facteurs p,: 


olivine ... p; 

augite . . .. p, « labrador 
riebeckite .. p, - andésine 
biotite . . .. p, - albite 
muscovite . . Pp, - orthose 


pot: quartz 


et que linfluence de la granularité et des micro- 
discontinuités est mésurée par le facteur g,: 


roche à grain grossier, . . . . « . 85 
moven, . . . - * . 8 
» » »» fin, “ . “ “ . . . . 80 


» » » 


Si on admet que les essais d'altération expé- 
rimentale nous fourniront les facteurs p, et g,, 
on obtiendra indice d'altérabilité à partir de la 
composition modale micropétrographique qui 
nous permet de calculer 4, si Von fait: 

a) le calcul de la composition modale de la 
roche, 

b) la construction du Tableau I qui va nous 
permettre de determiner k',.,, (1” indice d'alté- 
rabilité parciel ne se rapportant qu'à l'influence 
de la composition minéralogique). 

L'indice d'altérabilité parciel k',.. est calculé 
par: 

E: 
min. A+B 


£ 


o» (5) 
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Além de construtores de alternadores 
de potência superior a 2000 kVA e de motores 
de c.c.ou c.a. de potência superior a 1000 CY, 
possuimos equipamento especializado 
para a reparação de máquinas eléctricas rotativas 
de grande potência incluindo 
equilibragem dinâmica, ensaios de rigidez de 
bobines, perdas dieléctricas de isolantes, 
ring test e medidas de rendimento. 
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TABLEAU I 


Evolution de la roche dans 


Composition pondérée par 


CT ———— — e — 2 —— —— ———————e ee me em e 


Composition pondérée par 


sa oa nstettro p; des minéraux | les factures p; des minéraux Eai 
sains altérés 
ne 
te t 
o C é É ridd 

Roche à lorigine des temps e min =1 

t=t, 
Roche dans linstant géologi- A B 0< Kaio 1 
que dont on Vétudie | 

t= trinal 

Roche à la fin de la décompo- zéro Cc Kmin =0 


sition totale 


Note 1: On calcule B à partir des minéraux primaires qui ont donné les minéraux altéres, tout en les ponderant avec les 


convénients facteurs p;. 


Note2:C=C'=A+4B. 


Pour faire intervenir influence de la texture et 
des microdiscontinuités, voire de la surface éfec- 
tive de réactivité par unité de volume, on doit 
prendre les facteurs g,. 

Dans létat actuel de nos recherches on va 
admettre que: 


x=1,2,. « « « jusqu'ã 9 
On peut mettre que indice daltérabilité 


totale — k — será calculé par l'expression : 


== Kin — Bj mta Ei add (6) 
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Fourmaintraux (França), Skinner (Nova Zelândia), 
Struillou (França), Franklin (Inglaterra) e Guidicini 
(Brasil). 
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NOTICIÁRIO 


INTERKAMA 71 


Terá lugar, de 14 a 21 de Outubro de 1971 em 
Diisseldorf a Semana INTERKAMA — Congresso 
Internacional, com Exposição, da Técnica de Me- 
dida e de Automação. 

Os três pilares fundamentais da INTER- 
KAMA —71 são: 

1. Congresso científico-técnico 

2. Exposição especializada alargada 

3. Cursos práticos INTERKAMA, de partici- 
pação reduzida, consagrados à discussão, e a de- 
monstrações. 

Está desde já assente que o Congresso se 
ocupará sobretudo de problemas resultantes 
do emprego, cada vez mais generalizado, dos 
calculadores automáticos, cujas incidências são 
igualmente do domínio da técnica das medidas 
eléctricas de precisão. 

O prazo para inscrição dos expositores termina 
em Março de 1971. 

A 4º INTERKAMA realizada em 1968 teve 
702 expositores e 62000 visitantes, vindos de 
29 países, encontraram-se no terreno da exposi- 
ção. Quanto ao Congresso, registou 2500 parti- 
cipantes para as 56 conferências. Em relação aos 
cursos eles foram frequentados por 6000 pessoas. 


Conferência sobre Inteligência Artificial 


A conferência que terá lugar no Imperial College 
em Londres de 1 a 3 de Setembro de 1971, está 
a ser organizada pela British Computer Society 
em nome do International Joint Council on Arti- 
ficial Intelligence — uma organização formada na 
sequência da 1.º Conferência sobre Inteligência 
Artificial em Washington (U.S.A.) em 1968. 

Recebem-se comunicações sobre os seguintes 
temas : 

Fundamentos teóricos e estudos sobre inteli- 
gência artificial. 

Estudos teóricos de relêvo sobre ciências rela- 
cionadas com computadores, por exemplo, teoria 
matemática de programação. 

Demonstração de teoremas. 

Resolução de problemas heurísticos. 

Aprendizagem de máquinas. 

Pattern Recognition. 

Jogos. 
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“Compreensão” pelo computador por exemplo, 
comunicação em linguagem natural, por gráficos, 
fotografias ou outros meios de representação ; 
sistemas de pergunta-resposta. 

Pesquisas linguísticas relevantes para a inteli- 
gência artificial. 

Sistemas de inteligência artificial integrada 
por exemplo robots, sistemas de controle integrado, 

Hardware e Software relevantes para a inte- 
ligência artificial, por exemplo línguagens de alto 
nível, sistemas gráficos robot Hardware. 

Simbiose homem-máquina na resolução de 
problemas. 

Modelagem psíquica e fisiológica de sistemas 
inteligentes. 

Aplicações. 


Le No dd 
Aplicações do método PERT na indústria 
de construção 


“e mem o o mm e Amma, mm CS 


O INII vai realizar um curso sobre as Aplicações 
do Método PERT na Indústria de Construção. 

Este curso que terá lugar de 23 a 25 de Novem- 
bro de 1970 terá a duração de 18 h e será minis- 
trado em português pelo Eng.” Nelso Montes. 

O programa deste curso é o seguinte: 


— Actividades intervenientes nos empreendi- 
mentos de construção de edifícios. Programa 
de exigências funcionais. Análise técnico- 
-económica. Anteprojecto e projecto. Con- 
curso e adjudicação. Execução das obras. 

— Preparação da execução das obras. Esclare- 
cimento e harmonização dos projectos. 
Montagem do estaleiro. Elaboração dos 
programas de trabalhos, de pessoal e de 
máquinas. 

— Controle de execução da obra. Reuniões de 
coordenação. Controle dos programas. 

— PERT- Tempo. Princípios de base. Elaboração 
de listas. Definição de tempos. Estabele- 
cimento de redes. Cálculo automático de 
redes. Exemplo de aplicação. 

— PERT-Recursos. Princípios de base. Definição 
de recursos. Estabelecimento de cálculo 
automático de redes. Exemplo de aplicação. 


Horário : 9.30-12.30 e 14.30-17.00 h. 
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Tubos de caldeira, com e sem CORDÃO « do Levoxin 


Levoxin supera a protecção anticorrosiva com hidrazina. 
Considerada até hoje a mais eficaz. 


*Levoxin é um novo tipo de hidrazina 
com capacidade de reacção muito mais 
elevada, que proporciona uma segurança 
de serviço ainda maior nas suas insta- 
lações. 

Levoxin protege da corrosão todo o 
circuito de água-vapor mediante extrema 
aceleração da ligação de oxigénio, al- 
calinização, favorecimento da formação 
de camada protectora de magnetite, e 
graças ao efeito inibidor. 

É muito vasto o campo de aplicação 
deste produto Bayer: em caldeiras de 


toda a espécie, tubagens de vapor e 
de água condensada, aquecimentos a 
distância e por circulação, sistemas de 
refrigeração fechados. E também para a 
conservação, em húmido, de instalações 
que vão estar imobilizadas passageira- 
mente. 

Condições econômicas e segurança de 
serviço são os factores decisivos que 
falam a favor de Levoxin. Graças ao 
emprego de Levoxin, reduz-se a ocorrên- 
cia de grandes reparações e evita-se 
dispendiosas falhas de funcionamento. 


EESTI TT TETE TETE II ET ET TITO To 


Dicennsádio 778 AA 


: Queira enviar este cupão, colado sobre papel! 
: tUmbrado, à: t 
- BAYER PORTUGAL s.a.r.l. 

: Apartado 2365, LISBOA 

: Apartado 40, PORTO 


Pedimos o envio do prospecto sobre 
Levoxin. 


[| Pedimos conselhos para o seguinte 


problema: 
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TECNICA XVI 


GC. D. U. 621.373.4214 


AUTO-EXCITAÇÃO SIMULTÂNEA, ASSÍNCRONA E 

ESPONTÂNEA DE DOIS MODOS, NUM OSCILADOR 

DE REACÇÃO EMPREGANDO UM ÚNICO ELEMENTO 
NÃO-LINEAR ACTIVO DO 3.º GRAU 


RESUMO 


Faz-se a análise da cadeia fechada formada por dois 
ressoadores L OC ligados magnêticamente e um gerador de 
corrente comandado por tensão, cuja caracteristica é do 
tipo, 


= — au bu, a,5>0, 


O sistema de equações de evolução das amplitudes dos 
dois modos possíveis de oscilação, na hipótese assincrona, 
revela a existência dum estado de equilíbrio estável com 
amplitudes não-nulas de ambas as oscilações, que é aces- 
sivel a partir dum estado inicial de repouso. 

Estes resultados são confirmados experimentalmente por 
simulação analógica e ensaios num circuito real, 


I. INTRODUÇÃO 
A. Breve revisão dos resultados de estudos anteriores. 


O estudo da teoria não-linear da auto-excitação 
em sistemas com dois modos de oscilação pos- 
síveis, tem a sua origem em B. van der Pol [1] 
no início da década dos vinte. O sistema anali- 
sado por este Autor está esquematizado na Fig. 1 
e consiste num circuito passivo constituído por 
dois ressoadores ligados magnéticamente e um 
circuito activo não-linear formado pelo tríodo sob 
retroacção positiva em tensão, directa do circuito 
anódico para a grelha. O próprio Autor salienta 
que a teoria desenvolvida seria igualmente apli- 
cável a um sistema em que o elemento activo 
fosse do tipo «dinatrão», Fig. 2 (A). Grande parte 
dos estudos posteriores sobre este assunto têm 
incidido sobre o modelo da Fig. 2 (B) em que não 
figura ligação magnética. Com efeito, este por- 
menor não é essencial pois que no domínio das 


Por JF. BORGES DA SILVA 


Eng.” Electrotécnico Diplomado 
1.º Assistente LS, T. 


SYNOPSIS 


The closed loop formed by two magnetically coupled L O 
resonators and a voltage controlled current source charac- 
terised by a current-voltage relation, 


ab >0. 


i=—au+ bus, 


is analysed. The growth equations for the two possible 
oscullaticn mudes, assumed non-synchronous, show the exis- 
tence of a stable equilibrium with non-zero amplitudes for 
both modes, this equilibrium being attainable four a system 
initially at rest, The above results were confirmed by 
anatog aimulation and experiments involving a real 
circuit, 


oscilações quase-sinusoidais em sistemas quase- 
-Conservativos, apenas é relevante a existência de 
pares de polos conjugados do operador impe- 
dância do circuito passivo na vizinhança do eixo 
imaginário do plano das frequências complexas, 
independentemente da configuração particular do 
circuito utilizado. 

No caso presente, bastará pois que o circuito 
passivo apresente dois pares de polos conjugados, 
havendo ainda que exigir uma distância entre 
polos correspondentes dum e doutro par, sufi- 
cientemente grande comparada com a: distância 
desses polos ao eixo imaginário. Deste modo, com 
factores de qualidade dos ressoadores suficiente- 
mente elevados, o modelo da Fig. 2 (B) conduz 
aos mesmos resultados dos modelos anteriores. 

No seu estudo [1] B.van der Pol representa a 
característica não linear do elemento activo por 
um polinómio impar do 3.º grau, concluindo que, 
embora exista um estado de equilíbrio com am- 
plitudes não nulas dos dois modos considerados 


* Comunicação do Laboratório de Medidas Eléctricas do I.S.T. apresentada no Simpósio Sobre a Teoria da Infor- 
mação e dos Sistemas que teve lugar no 1.S.T. de 22 a 25 de Setembro 1970. 
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assincronos (frequências incomensuráveis), esse 
estado de equilíbrio é instável não sendo portanto 
realizável. 


Fig. 1 — Montagem referida por B.van der Pol no seu 
estudo /1) 


É geralmente atribuído a L. V. Skinner [5] ter 
constatado pela primeira vez a realizabilidade 
de oscilações simultâneas em duas frequências 


eee o 
EA 


START ET ET 
R.E. Fontana [6] concluindo-se da existência de 
estados estáveis, por exemplo, para uma relação 
2:1 das frequências mesmo com uma não lineri- 
dade do 3.º grau incluindo porém o termo quadrático. 

O progresso realizado no estudo destas matérias 
até 1954 pode ser apreciado num artigo de sintese 
devido a ]. S. Schaffner [7]. 

Mais recentemente, vários autores têm-se 
debruçado de novo sobre o problema das oscilações 
assincronas simultâneas em torno do modelo da 
Fig. 2 (B), corroborando resultados anteriores 
relativos às não-linearidades do 3.º e do 5.º grau, 
I. Bruyland [10], [13], H. ]J. Reich [11], [12]. 
Nestes artigos encontra-se também o estudo dum 
novo caso que consiste em ter, além do elemento 
activo não-linear do 3.º grau, dois elementos 
dissipativos não-lineares igualmente do 3.º grau, 
(i=au+by',a,b> o), cada um em paralelo 
com um dos ressoadores da Fig. 2 (B). O estado 
de equilíbrio estável de oscilação simultânea dos 
dois modos que assim se consegue é acessível a 
partir do repouso, isto é, é espontâneo. Referência 
a esta possibilidade tinha já sido feita por 
H. ]. Reich et al. [9]. 
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Fig. 2 — Modelos equivalentes para o estudo do problema das oscilações simul- 
tâneas. R, — Dípolo activo do tipo «dinatrão» 


incomensuráveis, sendo necessário para isso uma 
não-linearidade do 5.º grau com a forma conve- 
niente. Este estado de oscilação simultânea estável 
não é porém acessível a partir do repouso, depen- 
dendo de condições iniciais especiais. 

A oscilação simultânea em frequências comen- 
suráveis (caso sincrono), constitui um problema 
essencialmente diferente deste a que se refere o 
presente trabalho. Foi inicialmente estudado por 
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Outro caso ainda, conduz à mesma possibilidade 
e foi estudado por M. 1. Disman e W. A. Edson 
[8]. Consiste no emprego dum elemento activo 
com característica exponencial em conjunto com 
um mecanismo de variação automática do ponto 
de funcionamento em resposta à soma dos valores 
máximos das duas oscilações. Segundo os autores, 
o estudo é relevante em relação ao funcionamento 
dos osciladores de classe-C. 
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Finalmente e muito recentemente, certos pro- 
blemas desta índole, mas relativos a casos 
envolvendo um grande número de modos possíveis 
de oscilação, têm merecido também a atenção de 
alguns autores [14], [15]. 


B. Modelo do sistema investigado no presente estudo. 


Na Fig. 3 está representada a associação em 
cadeia fechada dum gerador de corrente ideal 
comandado por tensão e dois ressoadores quase- 
-conservativos ligados magneticamente. O circuito 
passivo receberá energia para a manutensão das 
oscilações através do par de terminais A, B, do 
ressoador 1, mas o gerador é comandado pela 
tensão aos terminais A, B, do ressoador 2. Este 
último facto torna o modelo da Fig. 3 essencial- 
mente diferente do da Fig. 1, não sendo já 
possível a sua redução aos modelos mais simples 
da Fig. 2. Com efeito, como é conhecido, os dois 
ressoadores oscilam praticamente em fase ou 
prâticamente em oposição de fase (rigorosamente 
no caso conservativo), consoante se considera a 
oscilação respectivamente do modo de frequência 
mais baixa ou mais alta. Assim a retroacção terá 
carácter positivo (desamortecimento) para um 
modo e carácter negativo (amortecimento) para 
o outro. Em qualquer dos modelos das Figs. 1 
e 2 0 comportamento do elemento activo é idêntico 
perante cada um dos modos considerados isolada- 


mente. 
Admitindo no modelo da Fig. 3 uma carac- 
terística i==f(u.) apresentando uma simples 


saturação, obtém-se o caso comum do oscilador 
de reacção por transformador duplamente sintoni- 
zado. 

Nestas circunstâncias, uma troca de ligações 
que inverta a tensão de comando «, em relação 


Fig. 3 — Modelo do sistema investigado 


TÉCNICA Nº 401 


à tensão », aos terminais do ressoador 1, permite 
obter ora um, ora outro dos dois modos possí- 
veis de oscilação do circuito passivo [3]. Um tal 
sistema nem sequer deve ser considerado como 
um oscilador com dois modos possíveis de osci- 
lação visto que, efectivamente, com uma dada 
configuração só um dos modos de oscilação é 
desamortecido. Note-se que este comportamento 
difere essencialmente do comportamento do 
modelo da Fig. 1 em que, embora não seja possi- 
vel obter simultâneamente as duas oscilações, 
é todavia possível sem mudança de configuração, 
obter uma ou outra por simples escolha de 
condições iniciais apropriadas. 

A possibilidade de oscilação simultânea em 
duas frequências surge no modelo de Fig. 3 se 
se admitir uma característica da forma, 


i=f(u)=—-au bu a,b>0 (1) 


representada na Fig. 4 (A). 

O elemento activo apresenta assim transmissão 
com inversão de fase no domínio das pequenas 
amplitudes, favorecendo a auto-excitação dum 
dos modos, e transmissão em fase no domínio 
das grandes amplitudes, favorecendo portanto 
o outro modo. A teoria e as experiências revelam 
então a existência de um único estado de equi- 
líbrio estável a que correspondem amplitudes de 
oscilação não nulas de ambos os modos. Esse 
estado é atingido espontâneamente a partir do 
repouso, começando por se auto-excitar o modo 
favorecido no domínio das pequenas amplitudes. 
O seu crescimento, levando o sistema ao domínio 
das grandes amplitudes, cria condições que 
possibilitam a auto-excitação do outro modo. 
O equilíbrio vem a ser atingido após sucessivas 


flutuações de amplitude alternadas dum e doutro 
modo. 

A forma da característica da Fig. 4 (A) também 
faz prever que se deva verificar o crescimento 
indefinido do modo favorecido pela transmissão 
em fase, desde que as condições iniciais sejam 
de molde a estabelecê-lo, logo de início, com 
amplitude suficientemente grande. 

Num sistema real é sempre de prever a exis- 
tência de saturação a partir de certa amplitude, 
como sugere a Fig. 4(B), a oscilação de grande 
amplitude mencionada acima virá pois, normal- 
mente, a ser estabilizada no domínio das grandes 
amplitudes. Nestas condições ter-se-á que uma 
perturbação inicial de grande amplitude conduz 


O circuito passivo fica assim caracterizado, com 
desprezo das perdas por agora, pelos valores de, 


L Cem, 


Encontrar-se-àá conveniência mais adiante em 
introduzir em vez destes os parâmetros, 


pk=——. 42) 


us G=y/€ M 


O circuito da Fig. 5 é então descrito, em ter- 
mos de valores instantâneos das tensões por, 


42 [1 k | | 
a -— tm = O 3 
dt? : | |: NES a 


Fig. 4 — Característica do elemento activo 


a um regime de grande amplitude com uma única 
frequência, aquela que corresponde ao modo 
favorecido pela transmissão em fase, enquanto 
que o arranque a partir do repouso ou provacado 
por uma perturbação de pequena amplitude 
conduzirá a um regime de oscilação simultânea 
nas duas frequências de amplitude moderada. 

Dado que a diferença de frequências dos dois 
modos é mormalmente pequena comparada com 
o próprio valor dessas frequências, estes regimes 
de oscilação simultânea nas duas frequências 
apresentam, no domínio do tempo, o aspecto 
característico de batimentos. 


II. EQUAÇÕES DE EVOLUÇÃO DAS AMEPLI- 
TUDES 


A. Equações do circuito passivo. 


Por simplicidade de análise, considera-se o caso 
em que os ressoadores são idênticos, Fig. 5. 
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Designando por u, e u, funções do tempo 
da forma, 


us=U,sen(r,t+ 2,) 


up == U, sen (vpt + Pp) 


em que , e «), são as raízes da equação caracte- 
rística do sistema (3), 


2 2 : 
( —  f, [o E. 
13 uy, k ' 


det | | |. |=0 


km” (3 migo 


Ly=Ly=L 


Cy=Co=C 


Fig. 5 — Circuito passivo 
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e que são dadas por, 


(1) n e a 


1oa= = do = aaa era f (5) 


VI +k 


a solução geral de (3) será da forma, 


idas (6) 


e a Rad O 


As soluções particulares linearmente indepen- 
dentes que compõem esta solução geral, associa- 
das cada uma a cada uma das raízes (5), 


aj Pl ad 1 ps 
Ho HA. Hs 2.5 


representam os dois modos de oscilação possíveis 
do circuito passivo. No seguimento, supor-se-á 
sempre serem as suas frequências incomensuráveis. 

O valor médio da energia eléctrica armazenada 
nos condensadores e correspondente a um dado 
modo de oscilação é dado por, 


A, = (,), + (A), = Cc — 


onde LI representa o valor máximo da oscilação 
do modo considerado, que é o mesmo num 
e noutro ressoador. Como consequência do sis- 
tema de equações (3), este valor médio é igual 
ao da energia magnética associada ao mesmo 
modo, pelo que a energia total armazenada em 


média será dada por, 
W=cu 


Para uma modificação lenta da amplitude 
à escala do período das oscilações, será neces- 
sário pôr em jogo uma potência, 


W=C— WU) (8) 


Nestas mesmas oircunstâncias, a existência de 
perdas qualquer que seja a sua origem, pode ser 
levada em linha de conta no balanço energético 
do circuito, referindo a potência correspondente 
ao valor máximo da tensão que caracteriza a 
amplitude do modo em causa através duma con- 
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dutância equivalente G. Acrescentando o efeito 
das perdas à expressão (8) obtém-se, 


Cc (u)+cu'=p (9) 


2 


onde P representa o valor médio da potência 
relativo ao fluxo de energia cedida ao circuito, 
associado ao modo de oscilação que se considera. 


B. Equações do circuito activo. 


De acordo com o esquema da Fig. 3, identifi- 
ca-se u, com —u,, resultando de (4) e (6) ser, 


U=—uo, =—u,stup=—U sento t+o,)+ 


+ Ussen (vost+9p) 


Utilizando a expressão (1) para caracterizar o 
circuito activo, procurem-se os termos da cor- 
rente i que apresentam as mesmas frequências 
que figuram na tensão u,. Encontra-se para os 
valores máximos respectivos expressões da forma, 


sa 13 2 . 
Lama Uq-3b (44 BA) 
4 2 


Ip=aUptab(DE q Sado) 
4 2 


(10) 


Os termos da corrente a que se referem estes 
valores máximos estão em fase (ou oposição) 
com os correspondentes termos da tensão u.. 
Na caracterização do elemento não-linear é 
mais informativo considerar, em vez dos coefi- 
cientes a e b, a inclinação na origem expressa 
através duma condutância G, e a abcissa das 
estacionaridades +-U,, Fig. 4 (A). Tem-se então, 


As parcelas da potência cedida ao circuito pas- 
sivo relativas a cada um dos modos serão então 
dadas por, 


(11) 


No valor médio da potência não figuram termos 
cruzados relativos a tensões e correntes de índice 
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diferente, por estes índices dizerem respeito a 
termos de frequências distintas. 


C. Equações de evolução normalizadas. 


Substituindo em expressões com a forma de 
(9), uma para cada modo, as expressões da potên- 
cia de (11) calculadas usando (10) obtém-se, 


[ed up=|-6,+6, 
de 

d 

CÊ fitje=) GA 
| dt | 


Observe-se nestas equações que a interacção 
dos dois modos condicionada pela não linearidade, 
se dá únicamente através das suas amplitudes, 
não figurando nestas equações as desfasagens em 
relação à origem dos tempos duma e doutra osci- 
lação. Esta circunstância é característica do trata- 
mento destes problemas na hipótese assincrona 
e resulta directamente de se admitir a incomen- 
surabilidade das frequências. Sob esta hipótese 
as desfasagens são irrelevantes para o balanço 
energético do sistema. 

Procede-se agora à normalização das equações 
(12): 

Toma-se G, como norma das condutâncias, 
introduzindo as condutâncias numéricas, 


(13) 


Introduz-se uma variável tempo numérica, 


e? —s (4) a 


0 


VLC (14) 


Finalmente normalizam-se as tensões tomando 
LU. para norma. Para variáveis dependentes do 
sistema de equações de evolução tomam-se então 
os múmeros a e b proporcionais aos quadrados 
dos valores máximos, 


O sistema de equações (12) escreve-se então, 


à GI 
| ta = . - —g,+ (1 — O ms DE 
dt" 26, 4 2 | 
ja, É 
js (1 er DE 
o (16) 
J 4 I 9 E. 
5 (— Us + á um [ui 
a o (12) 
E (+ U, + o | U,. 


A forma das equações (16) sugere uma nova 
normalização do tempo que absorva o factor em 
evidência no 2.º membro das duas equações. Na 
hipótese feita de evolução lenta à escala do 
período este factor tem de ser muito pequeno em 
relação à unidade. Introduz-se então nova variá- 
vel tempo numérica, 


+ = a (17) 


Nesta nova unidade o período das oscilações 
é expresso por um número muito pequeno 
comparado com a unidade, enquanto que na uni- 
dade definida por (14) era da ordem de 2r, 
A nova unidade está portanto proporcionada ao 
ritmo de evolução das amplitudes. 

Para atingir a máxima simplicidade na forma 
final das equações mormalizadas combinam-se 
ainda numa única constante em cada equação, 
os efeitos do amortecimento natural do circuito 
passivo com os efeitos de amortecimento ou 
desamortecimento provenientes do termo linear 
da característica do elemento activo. Definem-se 
para isso as duas constantes, 


a=4(1-—8,) 


E=4(1+8,) (18) 


As equações de (16) escrevem-se finalmente 
na forma, 


== (Dapa+aba 


ab d 
d 


(19) 
(—b+b+2a)b 
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Esta forma das equações de evolução é carac- 
terística do problema de duas oscilações assín- 
cronas ligadas por uma não linearidade do 3.º grau 
e surge com pequenas alterações em todas as 
diferentes variantes deste problema. A caracte- 
rística que distingue fundamentalmente o caso 
presente do caso da Fig. 1 e seus equivalentes, 
reside na ausência do sinal negativo no 2.º membro 
da 2.º equação ao contrário do que sucede nos 
outros casos. Esta troca de sinal provém directa- 
mente da existência de inversão de fase para o 
modo B ma transmissão através do circuito 
passivo, que condiciona a mudança de sinal da 
retroacção. Nos outros casos esta mudança de 
fase não tem óbviamente qualquer consequência. 

A dedução das equações de evolução das ampli- 
tudes foi conduzida segundo o método clássico 
aplicável aos sistemas quase-conservativos e que 
é fundamentalmente devido a B. van der Pol [2]. 
Este método é recomendável pela sua simplicidade 
e clareza no que respeita ao significado físico 
das aproximações introduzidas. Estes mesmos 
problemas são porém susceptíveis de tratamento 
matemático mais rigoroso, empregando o método 
de Krylcff e Bogoliuboff [4]. É porém conhecido 
que este método conduz nestes problemas aos 
mesmos resultados que o método clássico. A apli- 
cação do referido método a questões deste género 
pode ser apreciada nos trabalhos das referên- 
cias [6], [7], e |8]). A sua aplicação neste caso 
conduziu também aos mesmos resultados. 


HI. EQUILÍBRIO E ESTABILIDADE 
A. Pontos de equilíbrio. 


Igualando a zero os segundos membros de (19) 
obtém-se o sistema de equações de equilíbrio, 


O E ra es .() 
(a+2b—- a)a=€ (20) 
(b+2a-—-B)b=o 
cujas soluções dão as coordenadas dos pontos 
de equilíbrio, 
aj == O A, = & A; = 0 
(1) , (2) , (3) 
b/=o0 b; = 0 b, = f 
d a x — 24, — 2b, 
df | b | 2 bo 
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(28 —a) 
(4) (21) 


E 
-) 


(22— 


Ss a Us | ma 


O sistema (20) pode ser interpretado gráfica- 
mente com o auxílio das rectas representativas 
dos dois parêntesis igualados a zero. Esta 
construção, usada por vários autores [11], [12] 
e [13], está representada na Fig. 6. Nela, os 
pontos de equilíbrio estão indicados pela letra E 
seguida dum índice correspondente à numeração 
das soluções adoptada acima (21). As rectas 
figuradas são portanto isóclinas das trajectórias 
integrais do sistema (19) correspondentes às 
inclinações infinita e nula. Se se atender aos 
sinais algébricos das derivadas nas diferentes 
regiões do plano delimitadas pelas rectas da Fig. 6, 
imediatamente se reconhece a existência duma 


Fig. 6— Interpretação gráfica das equações de equilíbrio 


tendência das trajectórias integrais do sistema 
diferencial (19) para se enrolarem em torno 
do ponto E, . 


B. Estabilidade local. 


Representando genericamente por a, e b, os 
valores de equilíbrio de a e b, o sistema (19) 
linearizado ma imediata vizinhança dum ponto 
de equilíbrio toma a forma, 

-—= 24 


0 q 


(22) 
— +29, + 2b, 


emo b o 4 


Designando por A a matriz deste sistema, a 
equação característica do mesmo pode escrever-se, 


P—SAp+ dtA=0. (8) (23) 


A condição de estabilidade é a existência de 
todas as raízes no semi-plano real-negativo, que 
se exprime pelas condições, 


1) SpA< 0, 2) detA DO. 


A primeira condição é sempre satisfeita pois 
que, fazendo intervir (18), se tem, 


A =a-B=—4 (gy + 8p) 


e 4, € g, São necessariamente números positivos 
desde que haja dissipação no circuito passivo. 

Na discussão da segunda condição restringem-se 
as considerações aos únicos casos de interesse 
mediante as seguintes duas hipóteses: 


A) O sistema é efectivamente capaz de se 
auto-excitar no domínio das pequenas ampli- 
tudes através do modo A. 

Então, «0 ou o que é o mesmo, g,< 1. 

B) O estado de equilíbrio E, efectivamente 
existe, isto é, o ponto de intersecção das 
rectas da Fig. 6 ocorre no 1.º quadrante. 
Então atendendo a (21). 4, deve ser, 
a>0e b,>0,ou seja 268 >« e 2245. 
A primeira desigualdade verifica-se sempre, 
a segunda exige, 


Bo 
2 


84< (24) 


Conclui-se portanto que a hipótese B) acarreta 
também g,<<1, isto é, trocando a polaridade das 
ligações, Fig 3, também deve ser possível a auto- 
-excitação do modo B no domínio das pequenas 
amplitudes. A expressão (24) acarreta também 
para q, um condicionalismo mais forte do que 
resultava simplesmente de admitir a hipótese A), 
pois que exige g,< 1/2. 

Calculando agora det A para os diferentes 


(*) Sp A, traço de A. 
det À, determinante de A. 


pontos de equilíbrio, podem tirar-se conclusões 
quanto à sua estabilidade : 


Ponto E, detA = — «ai < U logo instável. 
Ponto E,, detA= — 3ab <O0 logo instável. 
Ponto E,, detA = —3%a,< 0 logo instável. 


Ponto E, detA= 37 b >0 logo estável! 


O único ponto de equilíbrio estável é pois 
aquele que corresponde à existência simultânea 
de duas oscilações. 

Analisando mais em pormenor este último 
ponto de equilíbrio, considere-se a situação limite 
em que o circuito passivo é conservativo. Ter- 
-Se-á 4, =$4p,=0, então 


a == [) == 4 (25) 


e em consequência Sp A = 0. As raizes da equa- 
ção (23) são nestas condições imaginárias puras. 
O ponto de equilíbrio E é portanto um centro. 
As coordenadas respectivas são, a —b = 4/3, 


Ocorre portanto sobre a bissectriz do 1.º qua- 
drante, Fig. 7. 


Existindo dissipação é sempre, 
a,—b=b-a=4 (gs, +39)>0 


logo o ponto E, ocorre abaixo da bissectriz do 
quadrante. As raízes da equação (23) são em 
princípio complexas, o que permite caracterizar 
E, como um foco estavel, A rapidez com que as 
trajectórias convergem para o foco é medida 
pela parte real dessas raízes tomada em módulo, 
ou seja por 


I 
-|Sp A | =2 (ga Rn) 


Porém, quando o ponto E, se situar abaixo da 
recta de equação, 


b=(7—V48)u 


Fig. 7, o discriminante da equação (23) torna-se 
positivo e as raízes ficam ambas reais negativas. 
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O ponto E, terá então o carácter dum nodo está- 
pe 


ty 


k 


Vig. 7 — Localização do ponto E , 


C. Separatrizes. 


Nos pontos de equilíbrio instável E,, E, e E, 
verifica-se que as raízes da equação (23) são 
reais, sendo uma positiva e a outra negativa. 
Estes pontos são, por consequência, do tipo colo. 
O par de raízes relativo a um ponto não é mais 
que o par de valores próprios da matriz A cal- 
culada nesse ponto, ao qual corresponde um par 
de direcções próprias que não são mais que as 
direcções das linhas separatrizes que se encon- 
tram nesse ponto de equilíbrio. A direcção pró- 
pria correspondente ao valor próprio positivo 
representa a separatriz divergente e a direcção 
própria correspondente ao valor próprio nega- 
tivo a separatriz convergente. 

A análise dos três pontos instáveis conduz aos 
seguintes resultados : 


/ . P, ntiios e. 0 
Valores próprios 
| os fo E O 
E, 
iai ur edge 
irecções próprias 
Ç Prop (0,1) 
| ,=— &<0 
Valores próprios 
|h=24-6>0 
a [ (1,0) 


Direcções próprias . 


( 2: e lô 
— 0º á 
| | 
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= t—- 2 EL] O 


p=B>o0 


Valores próprios 


E; 


ce. 0 
A od tido suaunte: é o - 
Direcções próprias 2 
h x (o , 1) 


Verifica-se no caso conservativo, x ==, que 
a separatriz que parte de E, é a mesma que 
atinge E, , sendo simplesmente uma recta de incli- 
nação — 1. No caso geral a separatriz que parte 
de E, vai terminar no ponto de equilíbrio está- 
vel E,, enquanto que aquela que termina em 
E, vem de (+ co, 0), Fig. 8. 


E, 2 


Fig. 8 — Separatrizes no caso geral 


IV. INTEGRAÇÃO DAS EQUAÇÕES DE 
EVOLUÇÃO 


A, Caso conservativo. 


Quando o circuito passivo não tem dissipação 
(25), o campo de velocidades definido polo sis- 
tema de equações de evolução (19) fica solenoidal 
puro, isto é, 


ov, dv, ' A h 
ie ; = DA PEÇE |, 2 7 ==. 
ia 9h («—2a— 2b) + (— É + 2b+4+24) 


Torna-se então possível definir as trajectórias 
integrais do sistema a partir duma função de 
fluxo & (a,b) apropriada, fazendo 


O (a,b) =.” 


A função de fluxo obtém-se por integração da 
diferencial. 


db=»v b da—v , dh 


que, em virtude de (26), é exacta. Obtém-se, 


P(ab)=ab(—-4 + ab) (27) 
onde a constante de integração foi escolhida de 
modo a tornar a função nula na origem o que 
arrasta ser nula sobre os eixos coordenados e 
sobre a recta 


-—4Hat-b=o0. 


A função, Fig. 9, é negativa no interior do 
triângulo E,, E,, E, e positiva no exterior. O 
ponto E, é um mínimo da função e fica envol- 
vido pelas trajectórias que nesta região são 
fechadas. Esta circunstância confirma a conclusão 
anterior, (III,B.), de que, no caso conservativo, 
E, é um centro. A esta afirmação de carácter 
local pode agora acrescentar-se outra de carácter 
global: Qualquer regime iniciado no interior do 
triângulo permanece, na sua evolução no tempo, 
indefinidamente no interior dum domínio limi- 
tado contido nesse triângulo. Pelo contrário, 
qualquer regime iniciado no exterior do triângulo 


Fig. 9 — Trajectórias no caso conservativo 
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crescerá indefinidamente segundo a amplitude do 
modo B. No caso conservativo, o interior do 
triângulo é portanto o domínio de estabilidade 
do sistema no que respeita às condições iniciais 
da sua evolução no tempo. 


B. Caso geral. 


Na impossibilidade de calcular as soluções por 
um método analítico simples, recorreu-se ao cál- 
culo analógico. A mecanização do sistema (19) 
foi feita segundo o esquema da Fig. 10, tendo 
sido obtido, entre outros, o traçado relativo ao 
caso de existir dissipação ligeira representado na 
Fig. 11. 

Como fora previsto, Fig. 6 e pela análise da 
estabilidade local de III, B., as trajectórias enro- 
lam-se em torno do ponto de equilíbrio antes de 
o atingirem. A existência de dissipação no cir- 
cuito passivo é portanto essencial para que o 
equilíbrio seja efectivamente atingido, pois que 
doutro modo, Fig. 9, as trajectórias são fechadas. 

O traçado da Fig. 11 revela ainda que o limite 


Mecanização das equações: 
A=-—-ara+ta+ 2h 


|--A , 


| - E mestã, Ds fdcidcgs 
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da região de atracção do ponto É, é a linha 
separatriz que termina em E,. Todas as trajectó- 
rias iniciadas sob essa linha convergem para E,. 
A região de atracção estende-se assim por uma 
faixa infinita, progressivamente mais estreita, 
disposta ao longo do eixo das abcissas. Ver 
também Fig. 8. Esta circunstância significa 
que condições iniciais; mesmo de grande 


andamento da separatriz que termina em E,. 
Este andamento pode sempre obter-se para cada 
caso, e à falta doutro método, pela integração 
numérica do sistema (19) partindo da vizinhança 
do ponto E, e na direcção conhecida da separatriz 
obtida em III, C. Dado que a separatriz termina 
em E,, é necessário inverter o sinal dos acrés- 
cimos de 7. 


Fig. 1x — Traçado das trajectórias obtido pela mecanização das equações de 
evolução segundo o esquema da Fig. ro, 


amplitude, desde que correspondam práticamente 
só à presença do modo A (modo favorecido no 
domínio das pequenas amplitudes), conduzem 
ainda assim à oscilação simultânea nas duas 
frequências. 

A chave do problema da estabilidade global 
deste sistema está portanto no conhecimento do 
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No Quadro I apresenta-se o resultado dum 
cálculo típico através dos valores das coordenadas 
de alguns pontos da separatriz. Os primeiros 
valores do quadro revelam um andamento prática- 
mente rectilíneo à saída de E,, os últimos valores 
revelam uma aproximação rápida do eixo das 
abcissas, ultrapassada a abcissa de E,. 


QUADRO i 
1 b MT, b 
a | q E 

0.0187 4.208 3.558 1.230 
Us 127 1.165 5. 438 0,417 
9,500 3.812 To 259 UV. 175 
1.108 3.247 12.10 0.040] 
2. 195 2.313 65.109 0. 0001 

a = 3.846 

à = 4.286 


Pontos da separatriz limite do domínio de atracção de E, 


V. EXPERIÊNCIAS 


A. Simulação analógica. 


Para confirmação da aplicabilidade da teoria 
feita ao modelo da Fig. 3, simulou-se esse 
modelo directamente com a montagem da Fig. 12, 
que não é mais que a mecanização do sistema 
de equações diferenciais dos valores instantâneos 
das tensões u, e u, da Fig. 3. Trata-se dum 


sistema de duas equações de 2.º ordem. À carac- 
terística do elemento activo não-linear foi apro- 
ximada por uma linha quebrada, Fig. 12 que 
traduz o essencial quanto à forma dessa 
característica. A simulação da Fig. 12 permite, 
além da observação do comportamento das ten- 
sões dos ressoadores u, e u, isolar também 
tensões relativas aos modos A e B para obser- 
vação, realizando a soma e a diferença das 
tensões dos dois ressoadores (6). 

Os oscilogramas n.º 1 a 7 foram obtidos cóm 
a montagem da Fig. 12. O grupo 1,2 e 3 corres- 
ponde ao arranque a partir do repouso. Tal como 
foi descrito na introdução, arranca primeiro 
o modo A (Osc. 1), depois o modo B (Osc. 2), 
antes de se atingir o regime estacionário 
notam-se sucessivas oscilações alternadas da 
amplitude dum e doutro modo (ignorar o den- 
teado das envolventes). O Osc. 3 mostra a 
sobreposição dos dois anteriores e representa 
a tensão do ressoador 1, o crescimento inicial 
deve-se ao modo A, o subsequente aparecimento 
do modo B dá origem ao batimento característico 
do regime estacionário final. O Osc. 4 apre- 
senta esse regime final em pormenor numa 
escala de tempos dilatada. O grupo de oscilo- 


”s “ 
RFestoador 1 


E nao- 


Ressoador 2 


Fig. 12 — Simulação do modelo da Fig. 3 
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se & raro apresentarmos 
o preco mais baixo, 
porque razão vendemos 
tanto ôleo para 
grandes turbinas? 


Das 46 turbinas a vapor, acima de 400MW, insla- 
ladas entre 1966 e 1968, 24º/, seleccionaram 
MOBIL desde que entraram em funcionamento. 

Estas vendas não se baseiam só na qualidade 
dos produtos. A qualidade do serviço também 
conta. À equipa técnica MOBIL acompanha cada 
passo do projecto-fazendo verificações com os 
construtores da turbina e da estrutura; seleccio- 
nando os lubrificantes adequados para o equipa- 
mento auxiliar; recomendando produtos anti- 
ferrugem para proteger as peças; assegurando 
a entrega de lubrificantes na hora exacta e, final- 
mente, treinando o pessoal que trabalhará na lubri- 
ficação. 

Para os clientes tudo isto representa menos 
preocupações, perdas de tempo e despesas. 

Garantindo este tipo de serviço a MOSBIL não 
E pode proporcionar, sempre, os preços mais baixos. 
das ca No entanto, os clientes estão mais interessados 
nos resultados a longo prazo do que nos custos 
iniciais. 


ma Fase toolado da manutenção das 
Empresas. Drpois de iniciado, porém, o mais natural 


Mobil 


vendemos mais vendendo menos 
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A Mobil si 
mostra-l 


inumeras maneiras 


vendemos mais vendendo menos 


SO 
Por ano poupamos à 
industria vários milhares 
de escudos. 
Uma grande parte 
deste dinheiro 
sai das nossas algibeiras. 


vendemos mais vendendo menos 


Os Mobiltac 
ventam-se . 
tanto em serviço, 


que dificilmente 
poderemos 


vendé-los muito. Eai 


vendemos mais vendendo menos 


vendemos mais vendendo menos 


Este simples ensaio 
mostra como reduzir 
os custos 


de manutenção 
dos sistemas hidráulicos. 


vendemos mais vendendo menos 


A industria de 
embalagem 
aprecia-nos 

porque poupamos... 
poupamos 

o seu dinheiro! 


vendemos mais vendendo menos 


Se é raro apresentarmos 
o preço mais baixo, 
porque razão vendemos 
tanto óleo para 

grandes turbinas? 


vendemos mais vendendo menos 


gramas 5, 6 e 7 por seu turno, representa pela 
mesma ordem o caso do arranque com amplitude 
inicial do modo B. Como este modo é amortecido 
no domínio das pequenas amplitudes, começa 
por desaparecer ((sc. 6), entretanto a sua 
presença impede a auto-excitação do modo À que 


por se observar o desaparecimento da oscilação 
inicial, devida ao modo B, para reaparecer mais 
tarde agora devida ao modo A. O reaparecimento 
do modo B inicia finalmente o batimento corres- 
pondente ao regime final. 

Na simulação a frequência das oscilações era 


OSCILOGRAMAS 


EVOLUÇÃO PARTINDO 
DO REPOUSO 


Mil 
| 


só vem a crescer depois do desaparecimento do 
modo B (Osc. 5). Uma vez alcançado o domínio 
das grandes amplitudes, através do crescimento 
do modo A, o modo B reaparece, (Osc. 5). 
Assim na tensão do ressoador (Osc. 7), começa 
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EVOLUÇÃO COM AMPLITUDE 
INICIAL DO MODO B 


de cerca de 15 Hz e a duração dos oscilogramas 
35 segundos! As condições são portanto efecti- 
vamente de evolução das amplitudes lenta à 
escala do período das oscilações. O factor de 
ligação simulado era de 0,055. 


B. Circuilo real. 


Com bobinas e condensadores constituiu-se 
um circuito passivo equivalente ao do modelo 
da Fig. 3 que se associou em cadeia fechada com 
um circuito activo cuja característica é da forma 
dada na Fig. 4(B). Ver Fig. 13. Para este fim 
empregou-se um amplificador de tensão diferen- 
cial (amplificação 40x), com uma resistência adi- 
cional em série na saída. Por meio do par de 
díodos associados à entrada + do amplificador 
consegue-se realizar a inflexão da característica 
no domínio das pequenas amplitudes. A limitação 
no domínio das grandes amplitudes é dada pela 
própria saturação do amplificador. 

O Osc. 8 representa um funcionamento típico 
em regime estacionário, no domínio das pequenas 
amplitudes com os dois modos de oscilação pre- 
sentes. O arranque a partir do repouso consegue- 
-se mediante um curto-circuito momentâneo dos 
pontos A e B, Fig. 13, pois que pela ligação 


cia de assincronismo. Informação semelhante 
pode ser obtida sincronizando a base de tempo 
com a oscilação e observando a envolvente. Este 
segundo método é o único que é imediatamente 
aplicável ao circuito real. Constitui também uma 
boa prova de assincronismo constatar a possibi- 
lidade de variação contínua e regular da dife- 
rença das frequências das oscilações produzidas 
com o factor de ligação k dos ressoadores. 

Experiências feitas não revelaram qualquer 
tendência para o sincronismo na montagem 
de simulação e ligeira tendência, na vizinhança 
de certos valores de k, no circuito real. Esta 
tendência explica-se por a relação de frequências 
no circuito real se poder aproximar de razões 
de números inteiros pequenos (Osc.8). Não sendo 
fácil ma montagem real aproximar a situação 
quase-conservativa como na montagem de simu- 
lação, é-se forçado a trabalhar com modos mais 
separados na frequência. Comparem-se os osci- 
logramas 4 e 8. 


Resistencias em K9 
Capacidades em nF 
Inductanciae em mE 


Fig. 13 — Esquema do circuito real ensaiado 


inicial da tensão contínua de alimentação do 
amplificador se obtém invariávelmente o regime 
de grande amplitude com uma única oscilação. 

As frequências presentes situavam-se na pro- 
ximidade dos 6KHz. 


D. Prova de assincronismo. 


Sincronizando a base de tempo do oscilógrafo 
com um dos modos e observando o outro, 
a ausência duma imagem estacionária, mesmo 
usando um alvo de longa persistência  evidên- 
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VI. CONCLUSÕES 


O estudo do modelo da Fig. 3 combas e nas 
equações de evolução das amplitudes dos dois 
modos de oscilação possíveis, revela a existência, 
em condições apropriadas, dum único estado 
de equilíbrio estável que corresponde à oscilação 
simultânea nas duas frequências próprias do cir- 
cuito passivo, sendo esse estado de equilíbrio 
acessível a partir do repouso. 

Em termos do sistema de equações de evolu- 
ção, esse estado de equilíbrio é um foco. A rapi- 
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dez de convergência é função da dissipação exis- 
tente no circuito passivo que tem portanto carác- 
ter essencial na manifestação deste regime de fun- 
cionamento. 

No que toca à não linearidade, a característica 
essencial que deve apresentar consiste em pro- 
porcionar transmissão em fase no domínio das 
grandes amplitudes e transmissão em oposição 
de fase no domínio das pequenas amplitudes, 
(ou o inverso). Com efeito foi possível confirmar 
experimentalmente os resultados da teoria numa 
simulação em que a não-linearidade era tradu- 
zida por uma linha quebrada em vez do polinó- 
mio ímpar do 3.º grau que explicitamente se admi- 
tiu na teoria. 

É também possível pôr em evidência a exis- 
tência de regimes de funcionamento deste género 
numa montagem experimental empregando um 
circuito passivo real, constituído por bobines 
e condensadores, desde que se realize um cir- 
cuito activo com a não-linearidade conveniente. 
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NOTICIÁRIO 


46.º Encontro anual da Associação Técnica 
de Pavimentos a Asfalto 


Realiza-se em Oklahoma City nos dias 15-17 
de Fevereiro de 1971 o 46.º encontro anual da 
Associação Técnica do Pavimento a Asfalto. 
O Encontro tem por finalidade o progresso de 
tecnologia do pavimento a asfalto através das 
pesquizas e comunicação entre os membros da 
Associação. 

Apresentam-se em seguida os temas sujeitos a 
discussão distribuída por 5 sessões : 


Sessão 1. BASES AND SUBGRADES ; DESIGN 
AND PERFORMANCE. L. L. Kole presiding 


Milan Krukar e J. €C. Cook, «Comparison of 
Washington State University Test Track Expe- 
rimental Pavement Rings + 2 and = 3». 


James A. Carroll e B. A. Brakey, «Experimental 
Work, Design and Construction of Asphalt 
Bases and Membranes in Colorado and Wyo- 
ming». 

T. Allen Haliburton, «Performance of Asphaltic 
Pavement Components e Systems on Expansive 
Oklahoma Subgrades». 

Chester McDowell e A. W. Smith, «Design; 
Control, and Interpretation of Tests for Bitu- 
minous Hot Mix Black Base Mixtures». 


Sessão 1]. MIXTURE PROPERTIES & PAVE- 
MENT ANALYSES. C. L. Monismith presiding 


Roger Cragg e Peter S. Pell, «The Dynamic 
Stiffness of Bituminous Road Materials». 


Richard A. Burgess, O. Kopvillem, e F. D. Young, 
«Ste. Anne Test Road — Relationships Between 
Predicted Fracture Temperatures and Low-Tem- 
perature Field Performance». 

John A. Deacon, «Equivalent Passages of Aircraft 
with Respect to Fatigue Distress of Flexible 
Airfield Pavements». 

Bruce G. Frenzel, Philip F. Dickson e Jack 5. 
Corlew, “Computer Analysis for Modification 


of Base Environmental Conditions to Permit 
Cold Weather Paving». 


Sessão III. SYMPOSIUM. «ASPHALT PAVING 
FOR THE 70'5”. V. M. Marker presiding 
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Sessão IV. PAVEMENT SURFACE CHARAC- 
TERISTICS. Ward K. Parr presiding 

A. Please e D. R. Lamb, «Binder Properties and 
the Retention of Texture in Asphaltic Surfa- 
cings». 

H. J. Fromm e J. T. Corkill, «Surface Course 
Mixes to Resist Wear by Studded Tires». 

J. H. Dillard e G. W. Maupin, Jr., «Use of Sprin- 
kle Treatment to Provide Skid Resistant Pave- 
ments». 

Gharles R. Marek, «Voids in Compacted Surface 
Treatment Aggregate Layers». 


Sessão V. MATERIALS PROPERTIES. W. J. 
Halstead presiding 


W. Heukelom e P. W. O. Wijga, «Viscosity of 
Dispersions as Governed by Concentration and 
Rate of Shear». 

]. €. Robertson e J. R. Moore, «A Correlation of 
the Studies of Inverse GLC of Paving Asphalts 
with Paired Comparisons Performance Ana- 

lysis». 

R. A. Jiminez e J. D. Kriegh, «Setting and Dura- 
bility of Asphalt Emulsions». 

R. L. Terrele C. K. Wang, «Early Curing Beha- 
vior of Cement Modified Asphalt Emulsion 
Mixtures». 


Modificação das datas de realização da 
feira de Hannover 


Inaugurada anualmente no último sábado de 
Abril e durante 9 dias a feira de Hannover vê 
a partir de 1971 alterada as suas datas de reali- 
zação com o acordo dos expositores de forma a 
aumentar o número de dias úteis de exposição. 

A partir do próximo ano a inauguração será 
feita numa quarta-feira e o certame terminará na 
6.º feira da semana seguinte incluindo apenas 1 
fim de semana e fazendo-se o possível para que 
o dia 1 de Maio — Dia do Trabalho na Alema- 
nha — não seja incluído todos os anos no tempo 
útil da feira. 

Deste modo resumiram-se os dias inúteis da 
feira de 5 para 2 dias. 
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EROSÃO DE PRAIAS — PANORAMA GERAL 
SITUAÇÃO EM PORTUGAL 


RESUMO 


Define-se praia e apontam se os principais agentes natu- 
rais que sobre ela actuam, Refere-se a interdependência 
entre esses agentes e o regime aluvionar costeiro, Ápontam- 
-se 08 casos de equilibrio estático e dinâmico de praias, os 
casos em que uma dada situação de equilíbrio pode ser 
interrompida, e foca-se a consequência grave dessa quebra 
de equilibrio, que é a erosão costeira. 

Citam-se métodos usados na defesa contra a erosão 
costeira e conclui-se pela necessidade de caracterizar a 
fisiografia costeira para a escolha do método mais indicado 
em cada caso, Mostra-se a necessidade de proceder a obser- 
vações de campo para à caracterização da fisiografia cos- 
teira, referindo-se, a propósito, programas de observações 
em curso nos Festados Unidos da América e na Hulanda, 
Finalmente, comenta-se a situação em Portugal quanto à 
erosão de praias e justifica-se a necessidade de se formular 
e pôr em práteca um programa integral de observações de 
agentes e processos costeiros, para o que se sugere um 
esquema de procedimento, 


1— INTRODUÇÃO 
1.1 — Generalidades 


Por praia entende-se correntemente uma área 
de aluviões — areias ou calhaus — erguendo-se em 
talude mais ou menos suave acima de um plano 
de água. 

Convém no entanto, para melhor compreensão 
do que se vai dizer, generalizar a noção de 
praia. 

Num conceito mais geral, a praia é uma massa 
de aluviões, em contacto com uma massa de água, 
condicionada pelas flutuações que pode sofrer o 
nível médio da superfície livre daquela massa de 
água e pelas perturbações geradas na mesma 
superfície. A praia estender-se-á acima e abaixo 
da superfície do mar. O limite do lado de terra 
é normalmente caracterizado por uma mudança 
da natureza do material da praia. Por vezes, em 
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SYNOPSIS 


The term «beach» às defined and the principal natural 
agents influencing beaches are listed, The interdependence 
between these agents and the alluvial coastal regime is 
pointed out. The cases of static and dynamic equitibrium of 
beaches are quoted, as well as those in which a given situa- 
tion of equilibrium can be interrupted, stressing the severe 
consequences that may fellow this interruption, namely 
coastal erósion, 

The methods used in the protection against coastal ero- 
ston are mentioned and it is concluded that one has to 
charneterize the coastal phystography in order to choose the 
method best surted to each case, The need to carry out field 
observarions on the characteristres of coustal physiography 
às shown, some observation programs under way in the 
United States and in Holand being mentioned to the purpose, 

Finally the situation of beach erosion in Portugal is 
examined and the need to formulate and put into practice a 
thourough program of ecbservatron of coastal agents and 
processes às stressed, and a course of action suggested, 


praias marítimas, esse limite corresponde ao 
aparecimento de vegetação e coincide com a cota 
mais alta atingida pelo espraiamento das ondas. 

É conhecimento resultante da observação 
directa de todos os que, como banhistas, frequen- 
tam as praias que as aluviões que as constituem, 
em particular as areias, se encontram em cons- 
tante movimento, subindo e descendo a praia, 
havendo ainda zonas em que a água se apre- 
senta com quantidades apreciáveis de aluviões 
em suspensão. O constante movimento das alu- 
viões é uma característica das praias e nessa 
característica se baseia um critério de delimitação 
do trecho submarino da praia. Poderá, com base 
nesse critério, considerar-se a praia como termi- 
nando do lado do mar na zona onde os movi- 
mentos das aluviões, resultantes da acção das 
ondas, sejam julgados insignificantes. Ora, sabido 
quão variáveis são as características das ondas 
ao longo do tempo, fácil é reconhecer que o 


limite da praia é igualmente muito variável. 
A variabilidade desse limite será ainda função 
da variação do nível médio da água, em conse- 
quência das marés, por exemplo. No que se 
segue, vão considerar-se praias marítimas de 
areias sujeitas ou não a marés, e, como pertur- 
bações da superfície do mar, a ondulação de 
período médio (6 a 20s) gerada ao largo da praia, 
a que se sobrepõem ou não as vagas geradas 
pelo vento local, de período curto (1 a 65). Às 
perturbações da superfície do mar resultantes da 
presença isolada ou simultânea da ondulação do 
largo e das vagas locais dar-se-á a designação 
genérica de ondas. 

As ondas são os agentes principais da movi- 
mentação das areias. Estas movimentam-se quer 
em suspensão, quer por arrastamento junto ao 
fundo. O movimento em suspensão predomina 
num trecho entre a linha de rebentação das 
ondas e a linha em que o nível do mar encontra 
a praia, sendo maior a concentração de areias 
nas vizinhanças da rebentação, onde é maior a 
turbulência. 

Quanto às suas direcções, os movimentos clas- 
sificam-se em dois grupos extremos : movimentos 
transversais e movimentos longitudinais. 

Os movimentos transversais são aqueles que 
se geram quando as ondas atacam frontalmente 
a praia, ou seja, quando as cristas das ondas são 
paralelas às batimétricas. Nestas condições-limite, 
as areias movimentam-se em ambos os sentidos, 
numa direcção normal às batimétricas, com ou 
sem predominância em qualquer dos sentidos do 
movimento. 

Os movimentos longitudinais são os que se 
processam paralelamente às batimétricas, em 
qualquer dos sentidos, com ou sem resultante em 
qualquer deles. 

Tanto os movimentos transversais como os 
longitudinais são fundamentalmente ditados pelas 
características das ondas actuantes, quer directa- 
mente pelo movimento ondulatório, quer indirecta- 
mente pelas correntes resultantes do processo 
ondulatório. O movimento transversal puro 
pode-se afirmar que não existe na realidade, já 
que pressupõe condições de regularidade da praia 
e das ondas que a observação corrente desmente. 
Existe no entanto uma permanente tendência, 
nunca totalmente concretizada, para o seu estabele- 
cimento, no sentido em que as areias são 
movimentadas sempre segundo a direcção da 
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propagação das ondas, ou seja, segundo a normal 
às cristas. Como a praia reage aó ataque das 
ondas procurando adaptar-se às mesmas de modo 
a que as batimétricas se tornem paralelas às 
cristas das ondas, resulta a tendência para o 
estabelecimento dos movimentos transversais 
atrás definidos. A constante variabilidade das 
características das ondas — direcção de propagação, 
alturas e períodos — torna pràticamenteirrealizável 
a adaptação perfeita das batimétricas às cristas 
e, consequentemente, o estabelecimento de movi- 
mentos transversais puros. 


1.2 — Às ondas como agentes da movimentação 
das areias 


Foi dito que os movimentos das areias resultam 
quer directamente do movimento ondulatório em 
si, quer das correntes motivadas pelo processo 
ondulatório. Importa portanto fazer algumas 
considerações sobre o movimento ondulatório, 
para melhor compreensão dos movimentos aluvio- 
nares que dele resultam. 

As ondas que atacam uma praia são frequente- 
mente geradas a centenas ou milhares de quiló- 
metros de distância, em profundidades de centenas 
ou milhares de metros. Enquanto as ondas se 
propagam a profundidades suficientemente gran- 
des para que essa propagação não seja influenciada 
pelo fundo, diz-se que a propagação se faz a 
grandes profundidades. Considera-se que as 
profundidades deixam de ser grandes quando a 
profundidade relativa se torna inferior a 0,5, 
entendendo-se por profundidade relativa Oo 
quociente da profundidade d pelo comprimento L 
da onda. As ondas passam então a propagar-se 
em profundidades médias, diminuindo o seu 
comprimento L à medida que diminui a profun- 
didade. A passagem das ondas de grandes para 
médias profundidades é acompanhada de uma 
modificação das órbitas das partículas, de circulares 
a elípticas. Nas órbitas circulares os raios 
diminuem da superfície para o fundo, até se 


anularem à profundidade d ; - Nas órbitas 


elípticas, de eixo maior horizontal, os eixos 
diminuem da superfície para o fundo, onde o eixo 
vertical se anula. Junto ao fundo, portanto, as 
trajectórias das partículas são rectilineas. Por 
profundidades médias entendem-se aquelas em 
que d/L está compreendido entre 0,5 e 0,04. 
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A essas profundidades, a velocidade de propaga- 
ção das ondas é função de d e de L, pela teoria 
das ondas sinusoidais de pequena amplitude, 
enquanto a grandes profundidades ela depende 
só de L. Quando d/L < 0,04, entra-se na região 
das pequenas profundidades, onde a velocidade 
de propagação se torna independente de L, 
variando apenas com d. 

A velocidade das partículas do líquido no seu 
movimento orbital, na região de profundidades 
médias, é proporcional, para ondas do mesmo 
período, às dimensões dos eixos das elipses. 
Como estes diminuem da superfície para o fundo, 
as velocidades diminuem no mesmo sentido, sendo 
a diferença entre as velocidades à superfície e 
no fundo tanto menor quanto menor for a 
profundidade, já que os eixos da elipse variam 
exponencialmente com a profundidade. A velo- 
cidade máxima junto ao fundo será, portanto, 
função da velocidade máxima à superfície e 
função, para a mesma profundidade, da amplitude 
da onda à superfície. Compreende-se que, a 
partir de certo valor da velocidade junto ao fundo, 
as aluviões ali existentes possam entrar em 
movimento por arrastamento, descrevendo trajec- 
tórias rectilíneas. Se as órbitas das partículas 
fossem fechadas, ou, mais precisamente, se as 
órbitas fossem rigorosamente elípticas em torno 
de pontos fixos, as aluviões oscilariam em torno 
da sua inicial posição de repouso. Ora, a teoria 
prevê, e a experiência confirma-o, que as órbitas 
são ligeiramente abertas, isto é, que a mesma 
partícula tende a descrever a sua órbita em torno 
de pontos diferentes que avançam no sentido de 
propagação da onda. Por esta via, há uma 
tendência para arrastamento do material de fundo 
em direcção à terra e um transporte de massa 
de égua no mesmo sentido. O contínuo afluxo 
de água à praia por transporte de massa cria um 
desnível na superfície média da massa de água 
— mais alta junto à praia —, desnível esse que 
provoca uma corrente da terra para o mar que 
compensa o transporte de massa. Aquela corrente 
pode por sua vez dar-se com velocidade suficiente 
para promover o arrastamento das aluviões. 

O movimento ondulatório, inicialmente harmó- 
nico e estável, vai-se tornando menos estável 
à medida que diminuem as profundidades, até 
entrar em franca instabilidade que culmina com 
a rebentação das ondas a profundidades que 
variam com a altura e o período das ondas e 
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com o declive da praia. O processo de rebentação 
é muito turbulento, pondo em suspensão quanti- 
dades apreciáveis de areia. Essa areia em suspensão 
é agora mais facilmente transportável por cor- 
rentes que se estabeleçam no seio da massa 
líquida, entre elas a descrita corrente de compen- 
sação do transporte de massa. Mas além desta 
corrente de compensação, outras são geradas pela 
rebentação das ondas. 

Primeiramente, temos o chamado jacto de 
rebentação, que sobe a praia até se transformar 
em potencial a sua energia cinética, após o que 
se inicia a descida. Além desta corrente, outra 
aparece quando as ondas atacam obliquamente 
a praia. Nestas circunstâncias, gera-se uma 
corrente paralela à costa, concentrada entre a 
linha de costa * e a zona de rebentação e mantida 
pela componente da energia das ondas paralela 
à costa. Esta corrente, que na literatura anglo- 
-saxônica é conhecida por «longshore-current», 
é designada em certas regiões costeiras de Portugal 
por corredoura, designação que se preconiza. 
Estas duas correntes, jacto de rebentação e 
corredoura, são usualmente os principais agentes 
transportadores de areia, combinando-se a última, 
muitas vezes, com a corrente de compensação 
do transporte de massa. 

Na Fig. 1 esquematiza-se o caminhamento 
de areias resultantes das duas correntes. Acima 
da linha de costa, predomina o transporte 
em ziguezague resultante do jacto de rebentação. 


( Jacto de 
«TRANSPORTE EM ZIGUE ZAGUE rebentação). 
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7 TRANSPORTE LONGITUDINAL (Corredoura) 


Fig. 1 — Esquematização do transporte de areias 
pelas ondas 


No caso do ataque frontal da praia, não exis- 
tirá a corredoura e a movimentação das areias 
será condicionada pela infiuência relativa do trans- 
porte de massa e da sua corrente de compensa- 


* Linha de intersecção do nível médio da água em 
repouso com a costa, 
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ção, influência variável com as características 
da ondulação. Conforme predomina a primeira 
ou a segunda, a praia engorda ou é erodida. 

No caso do ataque oblíquo, a velocidade média 
da corredoura será função, fundamentalmente, 
do ângulo que a onda ao rebentar faz com a linha 
de costa, da altura e período da onda e da incli- 
nação da praia. 

A quantidade de areia de dada granulometria 
transportada por arrastamento ou em suspensão 
pela corredoura aumenta com a velocidade desta. 


2 — EVOLUÇÃO MORFOLÓGICA DE PRAIAS. 
FROSÃO DE PRAIAS 


2.1 — Generalidades 


Viu-se que uma característica das praias 
é a constante movimentação das aluviões — em 
particular areias — que as formam, em conse- 
quência da actuação de forças naturais, com espe- 
cial realce para aquelas provenientes da incidên- 
cia das ondas e que se traduzem na geração 
de correntes capazes de transportar, por arrasta- 
mento ou em suspensão, volumes consideráveis 
de areias. O sentido do transporte pode ser trans- 
versal — aproximando-se da normal à linha 
de costa—ou longitudinal, ao longo da mesma 
linha. O transporte transversal é um reflexo 
de correntes várias, entre as quais a correspon- 
dente ao transporte de massa em direcção à costa 
e a resultante corrente de compensação, da costa 
para o largo. O transporte longitudinal é feito 
pela corredoura, corrente que aparece quando 
a praia é atacada obliquamente pelas ondas. Para 
o objectivo final do presente trabalho, não 
se considera necessário descrever correntes trans- 
versais e longitudinais, devidas ou não às ondas, 
que aparecem em circunstâncias diferentes das 
que foram consideradas. 

Importa fundamentalmente ter presente que 
correntes transversais e longitudinais, que ocor- 
rem em variadas circunstâncias, arrastam conti- 
nuamente de qualquer trecho de uma praia quan- 
tidades consideráveis de areias, maiores ou meno- 
res consoante a intensidade daquelas correntes, 
função, por sua vez, das características das ondas 
actuantes. 

Nesta base, vai-se procurar explicar esquema- 
ticamente alguns casos de comportamento de 
praias quanto à sua evolução morfológica. 
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Consideram-se dois tipos principais de praias: 


a) Praias com movimentos transversais impor- 
tantes e movimentos longitudinais insigni- 
ficantes; 

b) Praias com movimentos transversais e movi- 
mentos longitudinais importantes. 


2.2 — Praias com movimentos longitudinais 
desprezáveis em relação aos movimen- 
tos transversais 


Em praias deste tipo, a variabilidade mais ou 
menos cíclica, ao longo de um período de tempo 
suficientemente longo, em geral um ano, das 
características da ondulação, provoca variações 
morfológicas da praia igualmente cíclicas. A um 
periodo de tempo de actuação de ondas em que 
predomina o transporte de areia para a praia — 
predominância do transporte de massa sobre 
a correspondente corrente de compensação — 
sucede outro em que a situação se inverte, pas- 
sando o transporte pela corrente de compensa- 
ção a ser o predominante. Na primeira situação 
a praia engorda, na segunda é erodida, redu- 
zindo-se a largura da sua parte emersa. À longo 
prazo, a largura da praia flutua em torno de um 
valor médio. Como exemplo, cita-se o caso da praia 
de Copacabana onde a amplitude de variação da 
largura chega a ser, num ano, da ordem dos qua- 
renta metros, para uma largura média de sessenta 
metros. 

Praias nestas condições não acarretam proble- 
mas de maior e, atendendo à constância do valor 
médio da sua largura, diz-se que estão em equi- 
líbrio estático. 


2.3 — Praias com movimentos transversais e 
longitudinais importantes 


Admite-se que os movimentos transversais, só 
por si, tenderiam a manter a praia em equilíbrio 
estático. Há agora, no respeitante aos movimen- 
tos longitudinais, que considerar dois casos: 


a) Correntes longitudinais unidireccionais ; 


b) Correntes longitudinais bidireccionais. 


2.31 — Correntes longitudinais unidireccionais 


Neste caso, considerando um dado sector da 
praia, sai desse sector em cada instante um certo 
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